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ABSTRACT 



Time evolutions whose infinitesimal generator is a fractional time derivative 
arise generally in the long time limit. Such fractional time evolutions are con- 
sidered here for random walks. An exact relationship is established between 
the fractional master equation and a separable continuous time random walk 
of the Montroll- Weiss type. The waiting time density can be expressed using 
a generalized Mittag-Lemer function. The first moment of the waiting den- 
sity does not exist. 
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1 Fractional Time Evolution 



A series of recent investigations |], |[ || ||, |] has found that, in a suitable 
long time limit, the macroscopic time evolution T l of a physical observable 
X(t) is given as a convolution of the form 



T t a X(t ) = Jx(t -s)h c 



s\ ds 
t) T 



(1) 



where t > 0, < a < 1 and 



h a (x) = —Hll ( - 



(0,1) 
(0,1/a) 



(2) 



is defined through its Mellin transform 



H$ x 



(0,1) 
(0,1/a) 



x s 1 dx 



T{s/a) 



(3) 



For a = 1 



^i(x) = lim h a (x) = S(x — 1) 



(4) 



and hence eq. flU) reduces to 



T*X(t ) =X(t -t) 



(5) 



a simple translation. For fixed a the operators T*, t > form a one-parameter 
semigroup on a Banach space (e.g. the space of all continuous functions 
vanishing at infinity). 
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Much interest in the semigroups T l a derives from the fact that their in- 
finitesimal generators are proprotional to 



mm) = hghjW = _1_ pl')-^-')^ ( 6) 



+o s r(— Qt) 7 S' 







i.e. to fractional time derivatives of order a 0]. This may be seen from 
Ref. 0] (p. 302), and the fact that the kernels h a (s/t)/t in eq. ([[]) are stable 
probability densities (see also ||). 

Derivatives of noninteger order < a < 1 with respect to time are in 
this way found to furnish an important generalization for the infinitesimal 
generator of the macroscopic time evolution of physical observables. 

Given the generality of the results above it is natural to consider spe- 
cific examples by replacing the integer order time derivative with a frac- 
tional derivative. A particularly interesting example is given by diffusion 
and random walks [§, ||. The fractional generalization of the familiar master 
equation for random walks turns out to be intimately related to the theory 
of continuous time random walks with power-law tails in the waiting time 
density [fTofl . It is the purpose of this note to exhibit an exact relationship 
between a fractional master equation and continuous time random walks 

Let me conclude this introduction with the remark that only a special 
subclass of separable continuous time random walks with algebraically de- 
caying waiting time distributions corresponds exactly to a generalized master 
equation with Riemann-Liouville-type fractional time derivative. It would be 
interesting to specify precisely whether, and under which conditions, other 
waiting time densities with a power-law tail lead to the same result O, O. 



2 Fractional Master Equation 

Consider a random walk and let p(r, t) denote the probability density to find 
the walker at the position f 6 l d at time t if it was at the origin r = 



2 



at time t = 0. The positions r G M. d may be discrete or continuous. The 
conventional master equation for the random walk reads 

d 

_p(f, t) = 22 w(f- rXr", t) (7) 



with initial condition p(f, 0) = 5fo. To establish the fractional generalization 
of this initial value problem one must respect the nonlocal nature of the 
fractional derivatives. This leads to the fractional master equation 

d% + p(r,t) = fo — L^- +5>(f - r>(r<,t) (8) 

where < a < 1, t > 0, and the initial condition p(f, 0) = 5^ has been 
incorporated. Note that the fractional transition rates w(r) have units of 
(l/time) a . They obey the relation J2r w i^) = 0- Applying a fractional 
integral of order a to equation (H) yields an integral equation 

p(r, t) =8?o + -}- fit - tT' 1 E *>(?- *W, (9) 



reminiscent of the integral equation for continuous time random walks [13 



m m m o, p . 



3 Continuous Time Random Walks 



The basic integral equation for separable continous time random walks de- 
scribes a random walker in continous time without correlation between its 



spatial and temporal behaviour, and reads [16, ft 



p(f, t) = 6 n $(t) + I M - t')Y A(f - r>(r* ',t') dt 1 '. 
Jo ^ 



(10) 



3 



Here, as in (|9|), p(f,t) denotes the probability density to find the walker at 
the position r e IR d at time t if it started from the origin f= at time t = 0. 
A(f) is the probability for a displacement f in each single step, and i/j(t) is the 
waiting time distribution giving the probability density for the occurrence of 
a time interval t between two consecutive steps. The transition probabilities 
obey XV Mr) = 1- The function $(t) in eq. (|i~0]) is the survival probability 
at the initial position. It is related to the waiting time distribution through 



$m = i - C xp(t')dt'. (ii) 

Jo 

Note that in equation ( |T0D the system is prepared with the walker at position 
r = 0, and it develops from there according to i[)(t). 



4 Relation between equations (|9|) and (TO) 



The similarity between equations (||) and (10) suggests that the former is a 



special case of the latter. To show that this is true let 

POO 

V>(u) = &W(t)}(u) = / e- ut i){t)dt (12) 

Jo 

denote the Laplace transform of ^(i) and write 

A(g) = nmm = e e ^ A w (13) 



for the Fourier transform of A(f). Then the Fourier-Laplace transform p(q, u) 
of the solution to (|10D is given as |13|, |14|, |15], 



= 1 1 -V>(m) ^ 

u 1 - V>(w)A(ife) l-V>(u)A(g) 1 j 
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where is the Laplace transform of the survival probability. 



Similarly the fractional master equation (|]) can be solved in Fourier- 
Laplace space with the result 



P{q,u) = — Tzr (15) 

u a — w(q) 

where w(q) is the Fourier transform of the kernel w(r) in(|9]). Eliminating 
p(q, u) between ([14]) and fll5|) gives the result 

i^M = A(£^i = c (16) 
u a i/j(u) w(q) 

where C is a constant. The last equality holds because the left hand side of 
the first equality is <f-independent while the right hand side is independent 
of u. 

From fll6|) it is seen that the fractional master equation characterized by 
the kernel w(r) and the order a corresponds to a special class of space time 
decoupled continuous time random walks characterized by A(r) and ip{t). 
This correspondence is given precisely as 

Mu) = 1 —— (17) 

^ K ' 1 + Cu a K ' 

and 



X(q) = 1+Cw(q) 

with the same constant C appearing in 
the correspondence defines the waiting 
constant while the structure function is 
the transition rates. 



(18) 

both equations. Not unexpectedly 
time distribution uniquely up to a 
related to the Fourier transform of 
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To invert the Laplace transformation in ([T7|) and exhibit the form of the 
waiting time density if)(t) in the time domain one rewrites eq. (|P7D 



" " ! — (-1)* _ a(fc+ i 



^--TTW^o^ u (19) 



and inverts the series term by term. This yields the result 



j-a—l oo -r , ±01 \ k 

i>{t-a,C) = — V / (--) . (20) 



First one notes that for a = 1 the result reduces to 



iP(t;l,C) = ^exp(-t/C) (21) 



the familiar exponential form. For < a < 1 the result is recognized as 



il>(t;c t ,C) = —E a , a [--). (22) 



where the function 



« = £f(STfl (23) 



is the generalized Mittag-Leffler function (T^ . 



The asymptotic behaviour of ip(t) for t — > is readily obtained by noting 
that E a ^ a {Q) = 1. Hence ip(t) behaves as 



V>(*) oc r l+a (24) 
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for t — > 0. For < a < 1 the waiting time density becomes singular at the 
origin. 

The asymptotic behaviour for large waiting times t — > oo is obtained from 



the asymptotic series expansion for the Mittag-Lemer function [19 



N -n 

^ (zH -£r^r 0(|z| ~ N) (25) 



valid for |arg(— z)\ < (1 — (a/2))7r| and 2 — > 00. It follows from this that 
E aj3 (—x) oc x~ 2 for x — > 00 and hence that 



V>(t) oc r 1_a (26) 

for large t — » 00 and < a < 1. This result shows that the waiting time 
distribution has an algebraic tail as it is usually assumed in the theory of 



continuous time random walks [EG, 21, E2L E3L Oil 



In summary it has been shown that the master equation with a Riemann- 
Liouville fractional time derivative of order a corresponds exactly to a con- 
tinuous time random walk whose waiting time density is related to the gen- 
eralized Mittag-Leffler function and exhibits a power-law tail. It would be 
interesting to know the precise conditions under which the more general class 
of separable continuous time random walks obeying eq. (EBf) "approximates" 



11] this result. 
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